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ABSTRACT 

Measurements of the Hubble constant to distant galaxy clusters using the 
Sunyaev-Zeldovich effect are systematically low in comparison to values obtained 
by other means. These measurements usually assume a spherical isothermal 
(3 model for the intracluster medium (ICM). Formation processes and recent 
mergers guarantee that real clusters are neither spherical nor isothermal. We 
present the results of a statistical analysis of temperature bias in H determi- 
nations in a sample of 27 numerically simulated X-ray clusters drawn from a 
ACDM model at z=0.5 (the sample is online at http://sca.ncsa.uiuc.edu). We 
employ adaptive mesh refinement to provide high resolution (15.6 h~ 1 kpc) in 
cluster cores which dominate the X-ray and radio signals (sample images on 
line at http://cosmos.ucsd.edu/~wenlin/SZ/sz_cluster.html). The clusters pos- 
sess a variety of shapes and merger states which are computed self-consistently 
in a cosmological framework assuming adiabatic gas dynamics. We derive the 
value of H by computing the angular diameter distance to each cluster along 
three orthogonal lines of sight. Fitting synthetic X-ray and y-parameter maps 
to the standard isothermal beta model, we find a broad, skewed distribution in 
/ = H (S Z) / H (actual) with a mean, median, and standard deviation of 0.89, 
0.83 and 0.32 respectively, where H (SZ) is the value of H derived by using 
Sunyaev-Zeldovich effect method and Ho(actual) is the value used in the cos- 
mological simulation. We find that the clusters' declining temperature profiles 
systematically lower estimates of H by 10% to 20%. The declining temperature 
profile of our adiabatic system is consistent with the result including radiative 
cooling and supernovae feedback (Loken et al. 2002). We thereby introduce 
a non-isothermal (3 model as an improvement. Applying the non-isothermal (3 
model to the refined sample with well-fitted temperature profile, the value of f 
improves 9% relative to the actual value. The study of the morphology and the 
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clumping effects conclude that these two factors combine to induce scatter in f 
of ± 30 %. 



1. INTRODUCTION 

The Sunyaev-Zeldovich effect (SZE) is a spectral distortion of the cosmic microwave 
background (CMB) due to the inverse Compton scattering of CMB photons by high energy 
electrons in the hot intracluster medium (ICM) (Sunyaev & Zeldovich 1970,1972). At low 
frequencies (< 218 GHz) the distortion appears to be a decrement and thereby the peak of 
CMB spectrum (~ 160 GHz) intensity is reduced. This effect has been suggested as an im- 
portant tool for studying cosmology and astrophyics (Sunyaev & Zeldovich 1980; Birkinshaw 
1999 for more references). One important application of observing the SZE of galaxy clusters 
is the determination of the Hubble constant (H ). By combining SZE data with X-ray data 
from a cluster of galaxies, one can measure the distance to the cluster and thereby derive 
the value of H under a given cosmological model. 

A list of estimates of Hq based on this method are collected by Birkinshaw (1999). 
There are also some new results with better data and analysis (eg. Reese et al. 2000; 
Mason et al. 2001; Pointecouteau et al. 2002). The value of H from these findings are 
systematically lower than by other methods (Table. 1). This systematic discrepancy is due 
to the lack of information about the intracluster gasdynamics, which forces observers to 
assume that clusters are spherically symmetric, isothermal and in hydrostatic equilibrium. 
One well-known model based on these assumptions and frequently used by observers is the 
isothermal f3 model (Cavaliere & Fusca-Femiano 1976, 1978). However, drastic improvements 
of the X-ray observations in the past decade have revealed the deficiencies of this model. One 
particular important fact that contradicts the isothermal f3 model is the presence of clumping 
and substructures in clusters which characterize a cluster merger, (eg. Jones & Forman 1984; 
Gomez et al. 1997; Bilton et al. 1998; Neumann et al. 1999; van Dokkum et al. 1999). 
Moreover, non-isothermal temperature variations and decreasing temperature profiles have 
begun to be observed (Roettiger et al. 1995; Markevitch et al. 1997; Wang 1998; Honda et 
al. 1998). 

To investigate the consquences of the isothermal (3 model for H , we perform a statistical 
analysis on 27 rich clusters drawn from a fully hydrodynamic simulation with 15.6/h kpc 
resolution at cluster cores. We consider the influence of each assumption but focus on the 
temperature structures of clusters. In this paper, we describe the theory of the isothermal 
f3 model and the calculation of H in §2 and numerical methods of simulations in §3. We 
first derive the values of H by applying the isothermal (3 model in §4.1. Some statistics of 



- 3- 



the 13 parameter are discussed in §4.2; the non-isothermality and asphericity of the model 
are elucidated in §4.3 and §4.4. We then introduce a non-isothermal f3 model to improve the 
analysis method in §5. We also include the study of cluster morphology and the clumping 
effect in §6 and §7 and conclude in §8. 



2. ISOTHERMAL (3 MODEL : SZE AND H 

It is convenient to introduce a parametrized model to account the scattered photons 
in the intracluster gas. One popular model is the isothermal j3 model, which assumes that 
the electron temperature T e is constant and the electron number density has a spherically 
symmetric distribution: 

n e {r) = n e0 (l + — ) 2 (I) 

where r c is the core radius of the cluster, n e0 is the central electron number density, and (3 is 
a parameter close to |. (Cavaliere & Fusco-Femiano 1976, 1978). Over decades, this model 
has been much used in fitting the structures of X-ray clusters (Savazin 1988). If the electron 
density in a cluster atmosphere is n e (r), the Comptonization parameter along a particular 
line of sight is given by : 

/ ,k B T e 

cr T n e (r) -dl 

m e c z 

= D A [ a T n e {r)^di (2) 
J m e c 2 

where ks is the Boltzmann constant, ot is the Thompson cross section, m e is the electron 
mass, and c is the speed of light; the line of sight I is given by D A £, where Da is the angular 
diameter distance to the cluster and £ is the corresponding renormalized parameter. 

It can be shown that the X-ray surface brightness is given by : 
S x = An ^ +Z y J n e (r)n H (r)A eH dl 

= 4^i + z y DA J ne(r)n H (r)A eH d£, (3) 

where S x is the bolometric X-ray surface brightness, z is the redshift of the cluster, nn 
and n e are the hydrogen and electron number densities, and A e # is the X-ray cooling function. 
The maxima of the X-ray surface brightness and SZ y-parameter maps are denoted by S x q 
and z/q. By inserting Eqn.(l) to Eqn.(2) and Eqn.(3), we can eliminate n e o from Eqn.(2) and 
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Eqn. (3) and thereby solve for the angular diameter distance : 

= Vl W 2 ^W(3/3-0.5) £(|g) 1 

A " 47r(l + ^ l A; B T e j ^4 r(3/3) l r(|/3-0.5) J fl c 



where T(x) is the Gamma function and 8 C is the core size of the cluster. Since the angular 
diameter distance Da is a function of H and the cosmology, H can be calculated under a 
given cosmological model (Birkinshaw 1999; Reese 2000). 



3. NUMERICAL METHOD 

3.1. Cosmological Simulation 

The simulated clusters used in this project are part of our numerical "catalog" of clusters 
which are available at http://sca.ncsa.uiuc.edu (Simulated Cluster Archive; Norman et al. 
1999). The cosmological simulation is performed on the supercomputer Origin2000 at NCSA 
using our adaptive mesh refinement (AMR) cosmological hydrodynamic code ENZO. ENZO 
incorperates a Lagrangian particle-mesh (PM) algorithm to follow the collisionless dark mat- 
ter. The equations of gas dynamics are calculated by the piecewise parabolic method (PPM; 
Collela & Woodward, 1984) which yields high fidelity in the ICM temperature and density 
distributions. These numerical techniques have been described in more detail previously 
(Bryan et al. 1995; Bryan et al. 1999; Norman & Bryan, 1999). 

The cosmological parameters used in this simulation are: Q = 0.3, Qb = 0.026, Q\ = 
0.7, h = 0.7 and a 8 = 0.928. The high resolution adaptively refined simulations are performed 
after an initial low-resolution simulation is used to identify clusters. The simulation starts 
in an 128 3 base grid (level grid; L0), 256/h Mpc box at redshift 30. We then put two 
subgrids (LI & L2) covering each cluster to refine the Lagrangian volume there. Within the 
L2 grid, up to five more levels (L3 — > L7) of refinements are automatically introduced during 
the simulation. From LI, each grid is twice the resoultion of the previous one and therefore 
the best spatial resolution is 15.6/h kpc. (Loken et al. 2002) 



3.2. Analysis Method 

The simulated X-ray clusters are all at redshift 0.5. We observe these clusters along three 
lines of sight parallel to the edges of the box. Then we create synthetic maps of bolometric 
X-ray surface brightness and SZ y-parameter by projecting the two dimensional maps along 
the observational lines of sight. The projected surface brightness and the y parameter can 
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be written as: 



(5) 



(6) 



where 9 is the angular size of the cluster in the sky, and 6 C is the corresponding core size; 
S x o and yo are peaks in the projected maps. We derive both profiles by calculating the 
annular average of the projected surface brightness and y-parameter maps. Applying a 
three-parameter x 2 fitting on the profiles to Eqn.(5) and Eqn.(6), we can determine the 
central values (S x0 or y Q ) as well as 6 C and (5. The f3 value and 6 C used in Eqn.(4) are 
from fitting the X-ray surface brightness since the X-ray signal is much stronger than the 
SZ signal and dominate the final result. The X-ray cooling function A e # is calculated by a 
Raymond-Smith code with metalicity 0.3 solar. 



We analyze 27 numerically simulated clusters using the isothermal (3 model and thereby 
derive the value of H . Applying the analysis procedure mentioned in §3.2, we compute 
the values of H along three orthogonal lines of sight: x, y, and z. The resulting values of 
H are normalized by 70, which is the value adopted in the cosmological simulation. We 
denote the normalized Hubble constant by f and find a broad, skewed distribution in f with 
a mean and median of 0.897 and 0.829, respectively. The standard deviation is 0.32 if we 
parametrize the f distribution by a Gaussian (Table. 2). In order to identify the origin of 
this bias in Hq, we first test the assumption of "isothermality" . This is done by setting the 
temperature of the entire cluster equal to the emission weighted temperature and creating 
the corresponding bolometric X-ray and the SZE y-parameter maps. The results of analyzing 
these 27 isothermal clusters are displayed in Table. 3, where the mean and median are 1.111 
and 1.093 while the standard deviation is 0.297. Comparing the results from these two groups 
of clusters, we conclude that the applying an "isothermal" assumption on a non-isothermal 
cluster lowers the value of H by around 20 % compared to a real isothermal cluster. It 
is also noticeable that the standard deviation of the f distribution is close in both cases, 
which implies that the broad scattering in f is not mainly caused by the non-isothermality 
of clusters. 



4. RESULTS OF THE ISOTHERMAL (3 MODEL 



4.1. Bias Due to the Isothermal Assumption 
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4.2. (3 Value Distribution 

Based on fitting the X-ray surface brightness maps, some previous simulations have 
shown that the mean value of (3 is around 0.8 (Metzler & Evrard, 1994) while the observa- 
tional values are approximately 0.6 (eg. Mauskope et. al., 2000). The mean value of (3 in 
our simulation are 0.6210, 0.6054, and 0.6130 along 3 distinct lines of sight and is 0.61 if we 
combine the 3 line-of-sight results (Fig. 3; Table. 4). This indicates that our cosmological 
hydrodynamic code produces more realistic results than other codes (P3MSPH code, Met- 
zler & Evrard, 1994). Another interesting result here is that we actually find a correlation 
between f and (3 (Fig. 4), which deserves more investigation. We also fit the (3 parameter 
from the SZ contour map as well as from the X-ray surface brightness map. The (3 value 
derived from the SZ data is typically higher than the one derived from X-ray data. Fig. 5 
shows the relation of (3 values derived from the two different methods. 



4.3. Statistics of the Non-isothermality 

A rich cluster possesses a high potential well at the central region where most dark 
matter concentrates. At hydrostatic equalibrium, the temperature follows the gravitational 
potential. Since the potential decreases from the central core to the outer region, a de- 
creasing temperature profile is expected (§4.5). We define T iso as < ^° > to describe the 
"non-isothermality" of a cluster, where T and <T X > are the central and emission-weighted 
temperatures, respectively. Thus T iso equals 1 when a cluster is exactly isothermal. We 
denote the bias of the normalized Hubble constant f as df iso , which is defined as ^ s °|~^j or ^ , 
where f(iso) and f(ori) are the normalized Hubble constant averaged over three lines of sight 
for the isothermal and the original clusters, respectively (§4.1). The strong correlation (cor- 
relation coefficient equals 0.73) between these two variables confirms that the huge bias is 
due to the "isothermal" assumption (see Fig. 6). 



4.4. Elongation and H 

To study the elongation effect on H , we define a characteristic length for a cluster at 
certain direction as: 

where dl integrates through the peak of a projected X-ray surface brightness map along the 
observational line of sight. The results of three different lines of sight are denoted by l cx , 
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l C:V , and l c>z . Since the size of a cluster varies, we normalize l CjX , l C}V , and l C:Z by the value 
Imean = a/^U + ^c, y + ^c,z- n the °th er hand, for each single cluster, the normalized Hubble 
constant along three observational lines of sight f x , f y , and f z are renormalized by their 
mean value f mean — § (fx+fy+fz)- Because f m ean varies in different clusters due to different 
temperature structures (§4.3), we reduce corresponding effects by renormalizing f by f mea n 
. Fig. 7 shows the strong negtive correlation (correlation coefficient equals -0.74) between 
j^— and . * , which indicates that for a given aspherical cluster, the observed value of H 

i' mean Jmean 

is smaller along the longer axis. This result is consistent with the study by Roettiger et al. 
(1997). 



5. IMPROVED METHOD : NON-ISOTHERMAL (3 MODEL 

Our study of the 27 numerically simulated clusters reveals an universal temperature 
profile. The best fitting formula for the temperature profile is : 

T(r)=T (l + -)- a (8) 
r c 

where T is the central temperature of a cluster and r c is the core radius. At high redshift, 
such as 0.5 in our case, a lot of clusters are still dynamically young and non-relaxed, which 
results in the presence of some temperature clumps inside a cluster and thereby an irregular 
temperature profile. Here we disregard this type of clusters and concentrate on clusters 
without temperature clumps. Our refined sample includes 16 clusters with well-fitted tem- 
perature profiles. The average value of a drawn from these 16 clusters is 0.56. Fig. 8 shows 
a sample cluster fitted by this formula. 

Based on the dramatic improvement of X-ray data from Chandra and XMM, we antici- 
pate getting temperature profiles of clusters in the near future. However, if the temperature 
profile is not available, one can still use non-isothermal f3 model by taking the temperature 
profile derived from our simulation. If the temperature profile of a cluster follows Eqn.(8), 
the central temperature T and emission-weighted temperature <T X > are linked by : 

< T x >= Jo R (9) 

J e(r)4-Kr 2 dr 

where R is the radial size of cluster and normally chosen to be the virial radius, which is 
around 1 Mpc for rich clusters; e(r) oc n e (r) nnij) A(r), where the bolometric emissivity 
A(r) oc T(r) 1 / 2 , n e (r) and n#(r) are the number density profiles of electrons and the protons. 
Applying this temperature fitting formula to the final 16 clusters, we improve the mean value 
of H by 9 % (Table. 5). The lack of dumpings in the temperature profiles indicates the 
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relaxation of the clusters, which results in the refined sample having less scattering in the f 
distribution than the original ones. 

One might argue that this decrease in the temperature profile is a result of the insuffi- 
cient physics put in the simulation of the adiabatic clusters which do not have any cooling. 
However, some recent simulations including radiative cooling and supernovae feedback also 
show declining temperature profiles. These temperature profiles actually drop even more 
steeply than the adiabatic case (Loken et al. 2002) and therefore the improvement benefits 
from the non- isothermal j3 model would be even more for the cooling clusters. 



6. CLUSTER MORPHOLOGY AND H 

In practice, clusters are aspherical while the (3 model assumes that they are spherically 
symmetric. To describe the morphology of clusters in a more accurate way, we choose an 
elliptical model. An ellipsoid has three principal axes whose lengths we denote A, B and 
C from the longest one to the shortest one. The lengths of the three principal axes are 
propotional to the square root of the eigenvalues of the inertial tensor matrix. We examine 
the oblateness or prolateness of clusters by calculating their E-P values (Thomas et al. 1998): 

_ 1 B 2 (A 2 -C 2 ) 

= 2B 2 C 2 + A 2 C 2 + A 2 B 2 ^ 

_ 1 B 2 C 2 - 2A 2 C 2 + A 2 B 2 

= 2 B 2 C 2 + A 2 C 2 + A 2 B 2 
On the ellipticity-prolateness (E,P) plane, the strictly prolate clusters fall on the P = -E 
line while the strictly oblate clusters fall on the P = E line. Fig. 9 shows the distribution of 
the original 27 clusters on the (E,P) plane with the slope of a linear fit equaling to -0.733. 
This indicates a strong prolateness among clusters. The prolateness of clusters is believed 
to be due to the cluster forming processes when they form from filaments. Similarly, Fig. 10 
displays the result of the 16 clusters with well-fitted temperature profiles (see §5) and the 
slope from a linear fit of these profiles gives -0.585. Comparing the results from these two 
groups of clusters, the larger slope of the 16 clusters implies that they are less prolate than 
the original 27 cluters. Since relaxed clusters appear to be more spherical (less prolate), this 
result is consistent with the selection done by their temperature profiles. 
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7. CLUMPING EFFECT 

Accretion shocks and merging effects during the history of cluster formation result in 
the small-scale density fluctuations of the intracluster medium (ICM). Therefore, all cluster 
atmospheres are expected to be clumpy. The dumpiness can be described by the clumping 
factor C n , which is defined as <™ e > 2 ; where n e is the electron number density. Then the true 
angular diameter distance is (Birkinshaw 1999): 

DAitrue) = C n D A^estimated) (12) 

This means that without accounting for the small-scale clumping, the value of H based 
on Eqn.(4) is larger than the actual value. In our result, the values of C n range from 1.0 
to 1.5 with a mean of 1.25 and a median of 1.17 (Table. 6), which agrees with the result of 
Mathiesen & Evrard (1999). Applying the clumping correction lowers the mean and median 
values of f while the scattering of the f distribution decreases notably. Fig. 11 shows the f 
distribution of the original 27 clusters after the clumping correction; Table. 7 displays the 
results of the isothermal clusters before and after the clumping correction. 



8. CONCLUSION 

Utilizing the combined analysis of SZE and X-ray data with a standard isothermal (3 
model to determine the value of Hq results in the mean value of the normalized Hubble 
constant f being underestimated by 10 % and the median by 17 % compared to the actual 
value. There is a broad scattering in the f distribution and the standard deviation is around 
32 % of the actual value if we parametrize the f distribution by a Gaussian. By setting 
the temperature of a whole cluster to be the emission-weighted temperature, the mean and 
median of f become 1.111 and 1.093, which is 10% higher than the actual value. However, 
this overestimate is offset by the clumping correction (Table. 7). 

The discovery of an universal temperature profile motivates us to propose a non-isothermal 
(3 model to improve the analysis procedure. We refine our sample to 16 clusters which have 
well-fitted temperature profiles and apply the non-isothermal f3 model to analyze them. 
Based on these 16 relaxed clusters, the mean value of f is underestimated by around 21 % 
if the isothermal (3 model is applied and by around 12 % if the non- isothermal (3 model is 
applied, which improves the mean value of f by 9 %. By eliminating clusters with abnormal 
temperature profiles, we also eliminate non-relaxed clusters. The scattering in f distribution 
of these relaxed clusters is reduced to around 15 % of the actual value. 

We conclude that the broad scattering in the f distribution is due to the cluster morphol- 
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ogy and the clumping effect. All our simulated clusters are at redshift 0.5 and a lot of them 
are dynamically young. We have shown the strong prolateness in the cluster morphology 
distribution of the original 27 clusters as well as the 16 refined clusters. The asphericity of 
clusters and random observational angles between the clusters' axes and observational line 
of sight result in the broad scattering in the f distribution. Another reason for the broad 
scattering is the dumpiness within the ICM. We have shown that with the correction of the 
clumping factor, the scattering in the f distribution decreases while the mean and median 
drop as well. 

Therefore, in the real observation, we suggest to pick up relaxed clusters by their tem- 
perature profiles once the temperature structures can be measured. By applying a non- 
isothermal (3 model to these clusters, we can derive the value of H . Under this method, we 
anticipate ± 15% scattering around the mean value of H while the mean value is believed 
to be around 12% lower than the actual value. 

This work was supported by NSF grant AST-9803137 and NASA grant NAG5-7404. 
All simulations were carried out on the SGI/CRAY Origin2000 at the NCSA. 
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Fig. 1. — A 2D projection of our numerically simulated cluster. The color map is the X-ray 
surface brightness; the intensity drops from the central red core to the outer blue region. 
The contour map is the SZ y-parameter in units of 10 -4 . 
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Fig. 2. — The distributions of the normalized Hubble constant f derived from the original 27 
clusters using the isothermal (3 model. The top three figures display the results observing 
along the x, y, and z directions, respectively; the bottom figure shows the result combining 
the three lines of sight data. 
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Fig. 3. — The distribution of (5 values derived from fitting the X-ray surface brightness. The 
top three figures show the results along the x, y, and z observational lines of sight and the 
bottom one combines the three lines of sight data. 
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Fig. 4. — The correlation between the normalized Hubble constant f and the value (3. The 
top three figures display the results of the original 27 clusters along x, y, and z directions, 
respectively. The bottom one shows the result by combining the three lines of sight data. 
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Fig. 5. — /3(X-ray) is the value of /3 derived from fitting the X-ray surface brightness. /3(SZ) 
is the value of (5 derived from the SZE y-parameter contour map. 
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Fig. 6. — Bias of f due to the non-isothermality. df iso is defined as ^ 3 °|~^j or '' > , where f(iso) is 
the result of the 27 isothermal clusters while f(ori) is the result from the original ones. Ti SO 
denotes the degree of non-isothermality and is defined as £° . 
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Fig. 7. — For each cluster, 



J_ 



fn 



represents the bias of f in each observational direction 



relative to the f mea n which is the averaged value along the three observational lines of sight 
for that given cluster. Similarly, , c denotes the bias of the characteristic length along 

'mean 

one observational line of sight relative to l m ean which is the characteristic size of the given 
cluster. 
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The temperature profile of cl0022 
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Fig. 8. — The temperature profile of a sample cluster. The black-solid curve denotes the 
temperature profile of the cluster cl0022; the red-dotted curve is the fitting result for the pro- 
file. The radial distance is normalized by the core radius of the cluster and the temperature 
is normalized by the emission weighted temperature. 
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Fig. 9. — The morphology distribution of the original 27 clusters. The slope of a linear fit is 
-0.733. 



-21 - 



0.25 



EP plane (16 clusters ) 



0.15 



0.05 



-0.05 




oblate 



slope - -0.585 



-0.15 



prolate 



-0.25 



0.1 0.2 0.3 

E 



0.4 0.5 



Fig. 10. — The morphology distribution of the 16 clusters with well-fitted temperature 
profiles. The slope of a linear fit is -0.585. 
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f distribution with clumping correction 




Fig. 11. — Distributions of the normalized Hubble constant f with corrections of the clumping 
effect. The top three figures are the results of the original 27 clusters observing from x, y 
and z directions, respectively. The bottom one is the result by combining the data of the 
three lines of sight. 
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Table 1: The history of the measurements of H 



year 


value 


Method 


Reference 


1988 


85 ± 10 


Tully-Fisher relation 


Pierce & Tully 


1994 


73 ± 7 


Type II supernovae and 
the expanding photoshpere method 


Schmide et al. 


1994 


80 ± 17 


Cepheid variable stars in the Virgo Cluster 
with space-based measurement 


Freedman et al. 


1994 


87 ± 7 


Cepheid variable stars in the Virgo Cluster 
with ground-based measurement 


Pierce et al. 


1995 


67 ± 7 


Type la - based value 


Riess, Press, & Kirshner 


1996 


55 ± 3 


Type la supernovae 


Schaefer 
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Table 2: The mean, median and standard deviation of the normalized Hubble constant for 
the original 27 clusters. f x , f y and f z denote the results observing from the x, y and z 
directions, respectively. 



/ 


fx fy fz fall 


mean 

median 

sig 


0.876 0.939 0.876 0.897 
0.763 0.891 0.8 0.829 
0.393 0.274 0.287 0.32 



Table 3: The mean, median and standard deviation of the normalized Hubble constant for 
the 27 isothermal clusters. Here f x , f y and f z denote the results observing from the x, y and 
z directions, respectively. 



/ 


fx fy fz fall 


mean 

median 

sig 


1.051 1.178 1.104 1.111 
0.966 1.193 1.119 1.093 
0.303 0.278 0.306 0.297 



Table 4: The mean, median, and standard deviation of (3 value distribution. (3 X , (3 y , and (5 Z 
are the results from the x, y, and z observational lines of sight. f3 a u is the result combining 
the three lines of sight data. 



p 


Px Py Pz Pall 


mean 

median 

sig 


0.621 0.605 0.613 0.613 
0.613 0.597 0.599 0.599 
0.144 0.134 0.146 0.140 
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Table 5: Comparing the results of the 16 clusters with well-fitted temperature profiles using 
the isothermal (3 model and the non- isothermal j3 model. 



model 


isothermal j3 model 


non-isothermal /3 model 


/ 


/ 


/ 


mean 


0.804 


0.883 


median 


0.796 


0.884 


sig 


0.150 


0.152 



Table 6: The mean, median, and standard deviation of the clumping factor distribution for 
the original 27 simulated clusters. 





clumping factor 


mean 

median 

sig 


1.254 
1.168 
2.929 



Table 7: Comparing the results of the 27 isothermal clusters before and after clumping 
correction using isothermal j3 model. 



model 


before clumping correction 


after clumping correction 


/ 


/ 


/ 


mean 


1.111 


0.91 


median 


1.093 


0.93 


sig 


0.297 


0.228 



